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Abstract 

Here we show explicitly the relationship between the functions used in 
the atomic pair distribution function (PDF) method and those commonly 
used in small angle scattering (SAS) analyses. The origin of the sloping 
baseline, — 4nrpo, in PDFs of bulk materials is identified as originating 
from the SAS intensity that is neglected in PDF measurements. The 
non-linear baseline in nanoparticles has the same origin, and contains 
information about the shape and size of the nanoparticles. 

1 Introduction 

The atomic pair distribution function (PDF) analysis of x-ray and neutron pow- 
der diffraction is growing in popularity with the advent of nanoscience and nan- 
otechnology. The technique is more than 70 years old [TJ [H [3] and was originally 
applied almost exclusively to the study of glass and amorphous structures [H [5l 
[HI El [HI - However, the approach is proving powerful in solving structure on the 
nanoscale [3], where traditional crystallographic methods break down [TUJ. In 
particular, the study of the structure of discrete nanoparticles using the PDF 
method has recently become a focus [III [IU H31 [HI HSl HZl HB EE HS1 EI] ■ 
The convergence of this new need with the availability of powerful sources of 
high energy synchrotron x-rays and spallation neutrons and fast computing is 
greatly expanding the power and applicability of the method. 

The PDF, G(r), is defined both as a function of the real-space pair density, 
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p(r), and the reciprocal space scattering, F(Q) = Q[S(Q) — 1], as follows: 

G{r) = - f F(Q)smQrdQ (1) 
^ Jo 

and 

G{r) = 4nr(p(r) - p ), (2) 

where p is the number density of the material [22 } 123 ] IM l25l 126 ] 127 ] 128 ] 
[29l0[3(3[3Tl[32[33[a[3^ Equation [2] in this form 

works well for bulk materials, but the negatively sloping baseline, —Anrpo, is no 
longer valid when the PDF is calculated from finite-sized objects such as discrete 
nanoparticles [16] . Motivated by the need for a rigorous definition of the form 
of this baseline we rederive these equations here. We show the correct form of 
Eq. [2] in a number of cases of practical interest such as discrete nanoparticles 
and nanoparticle PDFs calculated from bulk models. The important distinction 
is provided by the small angle scattering (SAS) intensity and we explicitly relate 
the commonly used PDF functions with commonly used results from SAS. We 
also show that the widely used pair of definitions for the PDF above are actually 
incompatible with each other and that the definition Eq. [1] does not give rise to 
Eq.[2]but rather to R(r)/r, where R(r) is the radial distribution function. This 
work therefore resolves a long-standing ambiguity in the PDF literature. 

Compared to the PDF of a bulk sample [5T], the PDF of a nanoparticle 
is attenuated with increasing-r by a function that is related to the form of 
the nanoparticle [42]. For simple shapes, such as spheroids, spherical shells, 
rods and discs, this nanoparticle form factor can be computed analytically [43] 
144] [45] [40] and integral equations exist for more complex shapes [46]. This 
lends itself to a simple nanoparticle modeling procedure where the nanoparticle 
PDF is calculated from the PDF of a bulk phase analogue by multiplying by 
the assumed nanoparticle form factor [42] [47] [20] . The approach is successful 
for extracting precise quantitative structural information about the crystalline 
core of nanoparticles, including defects and size-dependent bond-lengths [20] . 
and this functionality has been incorporated in the latest version of the PDF 
modeling software PDFgui [18]. However, the method is not applicable when the 
nanoparticle structure has no bulk phase analogue. This is the case in general, 
for example in nanoparticles with surface modifications [12| and inhomogeneous 
compositions such as core-shell nanoparticles [49]. In these cases, models of 
discrete nanoparticles must be applied. As we discuss below, this results in 
an ambiguity about the precise form of the measured correlation function, and 
therefore how to calculate it. Currently, this is dealt with quite successfully in 
an ad hoc way, as for example in |16[ I50[ I51[ 137] . This paper presents a rigorous 
definition of the form of the baseline in terms of the nanoparticle form factor. 

In Section [2] we rederive the equations giving rise to the PDF to show the 
precise relationship between the measured correlation function in an x-ray or 
neutron total scattering experiment and the underlying model. In Section [3] 
we make explicit the link between the commonly used PDF and small angle 
scattering equations. This has implications for calculating PDFs from discrete 
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nanoparticle models for quantitative comparison with data, which are discussed 
in Section [4] In Section [5] we discuss the conditions under which the PDF can 
be calculated in real space. Section \E\ contains a brief summary. 



2 Derivation of the PDF equations 

To understand the precise relationship between the commonly used PDF equa- 
tions, nanoparticle structures, and small angle scattering we rederive the PDF 
equations from the beginning since subtle details of the derivation that are 
often overlooked have a significant impact on discussion presented here. Fur- 
thermore, the full derivation is not reproduced even in many textbooks on the 
subject [2 [24l [3l [39] and so these subtleties are not widely appreciated in the 
community. We start from the scattering amplitude from a set of i atoms at 
points fi in the kinematical limit: 

^) = Y J ^)^ i 

(3) 

i 

If the scattering from these atoms were totally incoherent the total intensity 
would be the sum of the intensities from each atom, 

line = ^2 

i 

= E/*(Q)/ l (Q) 

i 

= X)^a/:(0)/a(Q) W 
a 

= Nj2^f*(Q)UQ) 

OL 

= N{f), 

where the sum over a is now over the different species of atoms in the sample 
with N being the total number of atoms, N a being the number of atoms of type 
a and where the concentration of species a is c a = N a /N. Similarly, we can 
define the sample-averaged scattering power, (/) = ^ Q c a f a and 

« (5) 

a(3 
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The full coherent scattering intensity is given hy which is 



1 3 - (6) 

Here we have dropped the Q-depcndcncc of the atomic scattering factors to 
simplify the notation, but the fs are understood to retain their Q-dcpcndence. 
We can separate out the self-scattering, i = j, for which fij = 0: 

'■ E#/i ■ E';'< ,: '' 5 "' 

™ - (7) 

= N(f 2 )+Y,f;fc iQ ' r7j , 

where we have used Eq. 51 resulting in an expression for the discrete scattering 
intensity for i ^ j as 

h = Ic~ N(f) 

y,.r:.r,<'" : " ' . (8) 

We want an expression for the total scattering structure function, S(Q), 
which is defined as jfrkz ~ ^-rjff^ - ■ The second term in this definition is the 
Laue monotonic diffuse scattering that comes about because of the imperfect 
cancellation of intensity at the destructive interference condition when atomic 
sites are occupied by atoms of different scattering strength. It results in a mono- 
tonic incoherent background even in the case of perfectly coherent scattering. 

To get S(Q) from Eq.[7]we therefore must normalize by the total number of 
scatterers, N, 

|=(/ 2 )+^E/;/^- ( 9 ) 

Subtracting the normalized self scattering term to get 



and then normalizing by (f) 2 , wc obtain 

Ic (P) 1 



N(.f) 2 (f) 2 N(f)i 
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Thus, 



I c (f 2 ) 

S{Q) 1 = nW w 
Id 

N(f) 2 

= N(f\2 X! h 



(12) 



This expression yields precisely S(Q) — 1 in terms of scattering from atoms in 
our sample. 

For an isotropic sample, e.g., a powder of crystals or nanoparticles, we as- 
sume there to be a crystallite with every orientation with equal probability and 
we can take an orientational average. Place the Q along z so that we can ex- 
press Q ■ r~Zj — Qrij cos 6. Then the orientational averaging means that 9 takes 
all values with equal probability. The sample-averaged intensity for a pair of 
atoms will therefore be 

— \T ^ ff? d6e [ ^ cos e rl sin 9 
r^d^d^.sintf 



\-Kr\AQr, 



«W (13) 



e -iQn. 



2iQrij 
sm(Qnj) 



Using this in Eq. [5] gives the average coherent scattering intensity, as expressed 
originally by Debye [52] . From this we get the total scattering structure function 
for an isotropic sample, 

(i4, 

Thus, the reduced total scattering structure function, F(Q) = Q[S(Q) — 1], is 

sin(Qry) 



(15) 



It is convenient at this point to get rid of the independence of the x-ray 
form- factors. They are assumed to be isotropic so depend only on Q and not 
Q, which is a good approximation for scattering from core electrons especially. 
Write f(Q) = /(0)/(Q), where f(Q) has value 1 at Q = and contains the in- 
dependence of the form-factor and f(0) w Z , where Z is the atomic number that 
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scales the form factor. The Morningstar- Warren approximation [5] is that the 
Q-dependent part of the form factors can be well approximated by an average 
Q-dependence, f(Q) = N 1 Yla c «/a(Q)- hi this case the Q-dependence, 



f(Q) , comes out of the double sums in Eq. [T5] on the top and the bottom 
and cancels out. The fs that remain are Q-independent, and normally replaced 
by the atomic number (modified by any anomalous scattering factors). The 
same result holds for neutron scattering where the fs are replaced by coherent 
neutron scattering lengths, b. These have no Q-dependence and therefore the 
approximate method for removing the Q-dependence is not needed. 

Now we want to consider the inverse Fourier transform of F(Q). Because 
F(Q) is an even function, we use the sine-Fourier transform, 

f(r) = - I F(Q)sm(Qr)dQ. (16) 
n Jo 

We choose the — prefactor so that the direct sine transform has a prefactor of 
1. This is precisely the definition of the PDF in Eq. [TJ From this we get 

<i\ 2 r 1 V- ft sin (Q^) ■ /nun 



2 f*f- r°° 



(17) 



which, if we confine ourselves to the positive axis only, is 

/W = WVT/P ? fif* ^ r - r «)- (18) 

We can interpret f(r) in terms of the radial distribution function (RDF). 
The RDF, denoted R(r), is defined for an elemental system such that for an 
arbitrary atom i at the origin, Ri(r)dr gives the number of atoms in a shell of 
thickness dr at a distance r from that atom and the total RDF is the average 
of the partial RDFs over each atom taken at the origin. Thus, the integral of 
the RDF between two bounds gives the number of atomic pairs per atom with 
separation within those bounds. Equation [17] yields this behavior if we multiply 
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by r. For a solid with a atomic species we get 

= JfpEEP (19) 

= TfV2 X] C "^I! /j ' 
" ' a jeS 

where 5* is the set of atoms with distance from atom i greater than a and less 
than b. In the case of just one atomic species, this reduces to 



r b f 2 
/ H(r)dr=^^l 

J a J ■ r- c 



as required. Thus, 



ies (20) 

iVa, 



'< p ) = — (21) 
= 47rrp(r). 

The second expression comes from the relationship between the RDF and the 
pair density. The pair density is defined such that J drdifidOr 2 sm(0) p(r) = 
J R(r)dr, so that Airr 2 p{r) = R(r). Comparing Equations [TBI and |2"T1 we see that 
the definition Eq.[U]docs not yield G(r) (Eq.© and strictly | / °° F(Q) sinQrdQ 
R(r)/r 7^ G(r). However, we see below that in practice it is G(r) and not R(r) 
that is obtained experimentally in most cases. 
Finally, reordering Eq. [2TJ we find 

^ ^ 4irr 

= 4nr 2 N(f) 2 ' fi ^ ~ ^ 

In reality, I C {Q) is measured down to a minimum Q due to the experimental 
setup. This means that in general the forward scattering contributions are lost. 
We will consider the impact of this on the measured real-space function / (r) . 
We rewrite the expression for the experimental / (r) as 



2 f°° 

f(r; Q mm ) = - / F(Q) sin(Qr)dQ 

2 rQmin 

= iirrp(r) - - / F(Q) sm(Qr)dQ. 
n Jo 



(23) 



Of course, we have a finite Q m ax as well, but this will be disregarded during the 
following discussion, as the effects arc well understood [53] . 
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3 Low angle scattering intensity 



We will now consider a number of explicit examples to understand how the 
missing forward scattering affects the measured f(r; Q m in)- 

Consider a bulk material of uniform density and infinite extent. Since it has 
no internal structure there is perfect cancellation of all the discrete scattering 
intensity, Id, everywhere except at Q = since for all wavelengths it is always 
possible to find a pair of volume elements where the scattering is exactly out of 
phase with each other and cancels. The definition of Id in terms of a double sum 
over atoms is not entirely appropriate for this case (since there are no atoms!), 
but it gives an intuitive feeling about the behavior. 

I d (Q) = (N 2 -N)(f) 2 S(Q), (24) 

which is a delta-function spike at Q = sitting on a zero background. The 
integrated intensity in the delta-function scales like the volume of the sample 
squared. Strictly speaking the sample has to be infinite in extent, and therefore 
N = oo, to get a perfect delta-function. From Eqs. [T2] and [2U we see that 
5(0) -1=N-1&N since lim Q ^ sin Qr/Qr = 1. 

Now we do this more formally. Assume the sample has a uniform number 
density, p . Since the scattering length is defined per atom, the scattering 
amplitude of a volume element df at position r is i(j = /?o(/)e I< ^' r 'dr. For an 
infinite crystal, the intensity is given by 

UQ) = P 2 (f) 2 I [^< f - p) dfdP, (25) 



where the integrals are over all space, which gives a delta- function at Q = 0. 

If the material is finite in extent then the integrals are finite. To evaluate 
the integrals, we can define a shape function s(r) such that inside the shape 
s = 1 and outside the shape, s = 0. For such a material, 

h{Q) = Pl{f? I I s{r) S {r 7 y^-^drdr>. (26) 



Let us do a change of variables so that r" = f—r 1 , and dr" = df, in which case 
we have 

Ic(Q) = pl{f) 2 I I s{?)s{? + r")e i<3 ' r "dPdr 7/ 

J J (27) 

= plif) 2 [ dr"e i(? - r7 ' / s{r')s(P +r 7 ')dr'. 



The second integral is a self convolution, or autocorrelation function, of the 
shape function. Let us define 

7o(r) = y f s(P)s{P + f)dP, (28) 



where V = J s(r)dr is the volume defined by the shape function. This ■joir) is 
the characteristic function of the shape [54] , and has been called the nanoparticle 
form factor in the PDF literature [IT] [HI [20] . Defined as such, -fo(r) has the 
following properties: 

70 (0) = 1 

7o(r)dr = V ^ 
This definition, and a convenient dropping of the double-primes gives 

UQ) = pKffV f lQ {f)e^ F dr. (30) 



In analogy with the discrete case we want to convert this to S(Q) — 1 



ic _ <r 

WIT 1 (fY 



■, and using the fact that N = poV we get 



= PoJ To^e^df- 



(/ 2 ) 

(31) 



The second term, Kp-> is very small compared to the first term. It is order 
unity, where the first term scales as N = poV, and can safely be ignored in 
most cases. 

Now we want to take the orientational average. This must be done with 
care as, in general, the orientation of the nanoparticle shape and the underlying 
structure are correlated. For example, the morphology of the particles (plates 
or needles) depends on easy growth directions of the underlying structure. As 
discussed by Gilbert [40j . this means that the shape function and the internal 
structure of the particle are not, in general, separable and it is not correct to 
get the scattered intensity by convolving the reciprocal-space intensity with the 
Fourier transform of the characteristic function. Things arc greatly simplified 
in the case where the underlying structure, or the nanoparticle shape, or both, 
are isotropic, or approximately so. Then we can denote the angle-averaged 
characteristic function as 

70(r) - 7 ° (r) - /d0/d^sin(g) • (32) 
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Using this and Eq. [TSlwe get 



3(Q)-l = p dr d<j> d6>7o(r)e i( 3v 2 sm6l 
Jo Jo Jo 

POO p2lT PIT 

= Po dr d(j> / d07o(r)e i< 3-*V 2 sin0 
Jo Jo Jo 



2jv 



(33) 



sin (J.r 2 



po dr d</> / d^7o(r) — — — r sin 



o Jo Jo 

sin Qr 2 
Po / 7o (rj — — 47rr dr. 



Qr 



Since the particles have no preferred orientation in space, we have broken the 
average of the product in the first line into the product of the averages. This 
gives 

/•CO 

F(Q)= 47rp r 7o (r)sin(Qr)dr. (34) 
Jo 

Noting that this is the direct sine-Fourier transform, we take the inverse trans- 
form to get 

2 f°° 

f u (r) = - F(Q)sm(Qr)dQ 

t Jo (35) 

= 4irp Q rj (r), 

where the subscript u indicates that this result is for a solid of uniform density 
distribution, po- 

Next we consider a macroscopic crystal. The difference in F(Q) is at higher 
Q where, instead of complete cancellation of all the discrete intensity it appears 
at distinct reciprocal lattice points as sharp Bragg peaks. Importantly, in the 
region of Q below the first Bragg peak, the distinct scattering is zero except 
at very low-Q where small angle scattering region is reached. The small angle 
scattering intensity, I sas from the crystal is identical to that from the solid with 
uniform density: I? t as = I^ystai- The small and wide angle scattering regions 
are well separated in Q and I sas decays to zero before Q m in is reached in the 
crystal. Thus, 



2 r Qmin 

fsas{r) - - / F(Q)sin(Qr)dQ 

JO 

= fu(r) 



(36) 



and therefore 



- / F(Q) sin(Qr)dQ = 4 W7o (r). 

Ti" Jo (37) 
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We are now in a position to understand in detail the nature of the measured 
PDF f(r;Q min ). Substituting Eq.[37]into Eq.[23]we get 

f(r;Q m in) = 4:irrp(r) - 4m-rp Jo(r)- (38) 

This is similar to the definition of the PDF from Eq. [51 except that 70 (r) appears 
in the sloping baseline term. 

4 Calculating f(r; Q m i n ) from models 

We can now consider the calculation of measured PDFs using Eq.[35]in a number 
of interesting limits. 

4.1 Calculating in real-space for bulk crystals 

In the case of bulk crystals, the region of interest in the PDF is usually r -C D, 
D being the smallest dimension of the crystal. In this region, 70 (r) ~ 1. Thus, 

f(r;Q min ) = G(r) 

= ^r(p bu i k (r) - po), 

which is the familiar definition of G(r) in Eq. [5] The pair density function, 
pbuik{r), is calculated from a model with periodic boundary conditions [55j [56] , 
or from a box of atoms that is much larger in extent than the range of r of 
interest |57j . using Eq. 1221 The average number density po is given by the 
number of atoms per unit volume, which in the case of crystals is the number 
of atoms in the unit cell divided by the unit cell volume. 

Two approaches are typically taken to account for thermal and zero-point 
motion of atoms. One approach is to assume that these motions are well ap- 
proximated by a Gaussian probability distribution, in which case the delta- 
functions may be convoluted by Gaussians of finite width to represent the mo- 
tion. If the motion is anisotropic, Gaussian distributions with different widths 
in different directions may be used. Because the function being convoluted is a 
delta- function, from a practical perspective Eq. [39] is simply modified so that a 
Guassian of the appropriate width is added instead of a delta-function [SSI EH] ■ 
In models with many thousands of atoms, which is sometimes the case for re- 
verse Monte Carlo methods, the probability distributions can be built up from 
the static ensemble itself and no convolution is carried out. In this case no 
presumption of Gaussian dynamics is made. In practice the Guassian approx- 
imation works very well in most cases, and deviations from Gaussian behavior 
can be accounted for by introducing disorder in the models, which is a convenient 
way of separating the harmonic and non-harmonic contributions to the struc- 
ture. The effects of correlated dynamics are accounted for using an r-dependent 
Gaussian broadening [5£l [58l [60] . 

A convolution is also often carried out to account for the termination effects 
of the Fourier transform. For example, if the data are simply terminated at 
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Qmax, as is often the case when data are available to high-Q, the model PDF 
must be convoluted with a sine function [3]. From Eqs. [TBI and [39l it is clear 
that the function that must be so convoluted is f (r; Q m in)- This may be done 
with an integral directly in real-space, but it is often quicker to use a fast 
Fourier transform and do the convolution as a product in reciprocal space taking 
advantage of the convolution theorem. Note that, because of the sloping baseline 
in/(r; Q m in), the two convolutions (thermal motion and termination effects) are 
applied at different times during the PDF calculation. The thermal convolution 
is applied to pbuik(r) (Eq.(22) and then this is converted to G(r) (Eq.[5]), to which 
is applied the termination convolution. Modeling programs [561 148) account for 
this correctly, but when individual peaks are fit in real-space to extract peak 
positions and intensities, this is often overlooked, though it seems reasonable to 
ignore the termination effects in most cases. 

4.2 Calculating in real-space for nanoparticles modeled as 
attenuated bulk crystals 

In this case Pbuik{f) is determined using a model of a bulk structure as de- 
scribed in Section |4~T1 The pair density, p(r), in Eq. [3H]is the function for the 
nanoparticle, which is approximated as jo{r)pb u ik{ r ) [42) . Thus, 

fir; Qmin) = ^rj (r){pbuik(r) - po). (40) 

This approach has been implemented in the PDFgui modeling software [48] 
and used successfully on rather well ordered CdSe nanocrystals [2D] . The main 
shortcoming is that effects that cannot be incorporated in the average structure, 
such as surface relaxations or core-shell inhomogeneities, cannot be modeled. 

If there is a distribution of nanoparticle sizes and shapes, the characteristic 
function, 70 (r) can be replaced with an appropriately averaged characteristic 
function 

l(r) = J To(r; Ri,R 2 , ■ ■ >(i?i, i?2, . . .)dRxdR 2 .... (41) 

Here, p(Ri, R 2 , ■ ■ ■ ) is the normalized distribution of nanoparticle shapes param- 
eterized by Ri, i?2, . • .. For example, for spherical nanoparticles of radius i?, 
p(Ri, R 2 , • ■ •) = p(R), the distribution of nanoparticle radii. Finally, we replace 
Eq. [40] with 

f(r;Q m in) = ^rj(r)(p bu i k (r) - p ). (42) 

Great care should be taken to ensure that the result is unique when refining a 
number of nanoparticle morphology parameters beyond one or two. 

4.3 Calculating as the Fourier transform of the properly 
normalized Debye Function 

This approach has been successfully used by a number of authors [TJl [61] . The 
F(Q) function is evaluated using Eq. 1 151 and then Fourier transformed to obtain 
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the desired real-space function. To account for thermal and zero-point motion 
in reciprocal-space calculations, Eq. [15] is replaced with a version that includes 
Dcbyc- Waller effects, 

Here, cr^- is the correlated broadening factor for the atom pair, which is the 
real-space Gaussian width discussed in the case of bulk crystals. 

As we show here, for a quantitative comparison with measured data, care 
must be taken with the Fourier transform so that it is carried out over the 
same range of Q as the experiment. The main drawback of the reciprocal-space 
approach is that it can be very slow compared to direct real-space calculation due 
to the long-range extent of the signal from each pair. It is generally preferred for 
smaller systems. However, due to recent algorithmic advances, the calculation of 
the Debye equation for larger systems can be greatly accelerated under certain 
circumstances [6T] , 

Using this method, the termination effects coming from the Fourier trans- 
form are implicity included provided the Fourier transforms to obtain the model 
and data PDFs are terminated with the same Q m in and Q Tl 



4.4 Calculating in real-space from discrete nanoparticle 
models 

In this case, Eq. [38] is used directly, where p(r) is calculated from a finite model 
of the discrete nanoparticle using Eq. 1221 The difficulty arises in determining 
a correct form for the baseline — 4:Trrpojo(r). Up until now, the shape of the 
baseline has been approximated using expansions of ad hoc mathematical func- 
tions [TH [50] [5TJ [37] . This is successful at approximating the behavior of the 
baseline. However, in this work we derive the explicit form of the baseline shape 
in terms of the characteristic function of the nanoparticle, the autocorrelation 
function of the nanoparticle shape. This suggests a number of approaches to 
calculate the PDF baseline in a more physical way. 

If we have accurate small angle scattering data from the samples, from Eq. 1381 
we see that we can compute the PDF baseline from the measured SAS via a 
Fourier transform. However, care must be exercised as the derivation assumes 
that the sample is made up of discrete nanoparticles. In general, clusters and 
aggregates of nanoparticles will form and small angle scattering signals from 
these structures on different length-scales will be present and must be separated. 
Also, scattering density fluctuations of any sort in the sample will affect the 
SAS signal, as discussed elsewhere [55]. None of these effects need be explicitly 
considered if the SAS signal is not used in the PDF definition, as in Eq. |3"51 
though an alternative method is then required to determine the baseline. 

The inclusion of SAS data has the potential to add significant value to any 
refinement of nanoparticle models from the PDF. Both small and large angle 
scattering contain information about the shape and size of nanoparticles, but 
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this information is nearly decoupled from the internal nanoparticle structure in 
the small angle scattering. This same information is in the PDF, but it can be 
obscured by structure features, such as when the PDF prematurely attenuates 
due to a complex or amorphous surface structure. In this case the nanoparticle 
size obtained from SAS will be larger than the apparent NP size obtained from 
the PDF that reflects the size of the coherent core structure. Even without the 
inclusion of SAS data into PDF refinements, we can can learn much from SAS 
analysis techniques. PDF nanoparticle refinements usually start with a simple 
model that includes atomic positions restricted by a shape. Therefore, PDF 
analysis can benefit from the various ab initio methods |63l I64| for determining 
the shape of a scatterer from the SAS. 

Without the use of the small angle intensity for determining 70(f), we can 
consider approaches to determine it self-consistently from the model since it is 
the autocorrelation of the particle shape which is directly available from the 
model itself by determining a "shrink-wrapping" of the atomistic model. On 
the contrary, when the internal structure is well known, but the size and shape 
distribution of nanoparticles is not, then the characteristic function can be pa- 
rameterized and refined to obtain the approximate nanoparticle dimensions as 
was done in [2D] . 



5 The extent of small angle scattering 

Wc have considered the two asymptotic situations here of including or excluding 
all the SAS. If the SAS is retained in the intensity that is Fourier transformed, 
the resulting real-space function obtained is Anrp(r) = R(r)/r (Eq.[T7]). Exclud- 
ing it all results in G(r) (Eq. [38)) . We now consider the possibility that some, 
but not all, of the small angle scattering is included in the Fourier transform. 
This might occur in the case of very small nanoparticles, for example, when the 
SAS extends to wider angles. Here we estimate the circumstances under which 
a significant amount of small angle intensity will appear in a wide angle PDF 
experiment for the case of a sphere of uniform density. The scattering intensity 
is given by [43] 

/(0)CX (Qlje iMQR) - QRcos(QR)}\ (44) 

where R is the radius of the sphere. By integrating this equation from to 
Qmin, and dividing by the total integrated intensity, we get an expression for 
the proportion of small angle intensity above Q m i n for a given nanoparticle 
diameter [4"5] : 

l ( x ) = l-J-^[(2x i -x 2 + 3)cos(2x) 

27rr a (45) 

+ x(x 2 + 6) sin(2x) + 4x 5 Si(2x) - (5a; 2 + 3) ] . 

Here, x = Q m inR and Si(a;) represents the sine integral, Si{x) = J? sm jf ^ dx' . 
For a very small nanoparticle of radius 5A, we see that < 0.01, correspond- 
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ing to Qmin ~ lA , which is typical for a RAPDF [H5] experiment. Thus, 
for even quite small nanoparticles, practically all small angle intensity is below 
Qmin = lA -1 and Eq. [38] is appropriate. However, care should be taken not to 
extend Qmin too low in Q in a measurement of a nanoparticulate system. 

In the few-atom limit, such as the case of discrete small molecules, the small 
and wide angle scattering are not cleanly separated. To produce a complete 
real-space signal, one can approximate the small angle scattering from a candi- 
date structure model. This approach is commonly used in the study of small 
molecules in the gas phase [55] . The Fourier transform of the estimated scatter- 
ing approximates R(r)/r, the nominal "experimental" or "modified" RDF. An 
equivalent method for obtaining the modified RDF from wide angle scattering 
alone is to add a baseline estimated from a model structure to f(r, Q m i n ) [67] . 
This approach has been successful for calculating the R(r)/r for many-atom 
nanoparticles. 

6 Summary 

The PDF is a valuable tool for identifying the form and the interior compo- 
sition of nanoscale materials. Whereas the oscillating component of the PDF 
gives information about the interatomic distances within the material, the PDF 
baseline is a function of the characteristic function, a measure of nanoparticle 
shape which has its origin in the SAS that is usually disregarded in a powder 
diffraction experiment. This characteristic function goes unnoticed in macro- 
scopic particles, where the PDF is observed at distances that are much smaller 
than the particle diameter. For nanoparticles, the PDF baseline, and therefore 
the characteristic function, cannot be disregarded. We have presented a full 
derivation of the PDF equation taking into account the missing SAS and have 
reviewed different methods for calculating the PDF for nanoparticles. Given 
the relationship between the PDF and SAS equations, there is potential benefit 
in incorporating SAS data and analysis methods into PDF studies. 
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